arXiv: 1508.04075vl [cond-mat.mes-hall] 17 Aug 2015 
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We have studied a Superconducting Quantum Interference SQUID device made from a single 
layer thin film of superconducting silicon. The superconducting layer is obtained by heavily doping 
a silicon wafer with boron atoms using the Gas Immersion Laser Doping (GILD) technique. The 
SQUID device is composed of two nano-bridges (Dayem bridges) in a loop and shows magnetic 
flux modulation at low temperature and low magnetic field. The overall behavior shows very good 
agreement with numerical simulations based on the Ginzburg-Landau equations. 


Superconducting devices are currently used in many 
applications where they have shown their great poten¬ 
tial. This ranges from detectors for astrophysics up to 
a single photon level detection |pQ to very sensitive mag¬ 
netic sensors with Superconducting QUantum Interfer¬ 
ence Devices SQUIDs [2]. Another promising and rapidly 
growing field of research where superconducting devices 
play an important role is quantum information. So far, 
the most advanced circuits for quantum computing are 
based on superconducting qubits [5]. Many efforts have 
been made to increase the number of qubit operations 
that can be made within the quantum coherence time, 
including improvements using quantum error correction 
schemes and quantum non destructive measurements [ 3 ]. 
Even though multiple coupled qubit manipulations have 
already been performed SI E|, the next step is to demon¬ 
strate the possibility of performing operations that in¬ 
volve multiple logical qubits. Another candidate for 
quantum information is the spin qubit |B|. First real¬ 
izations were obtained with two dimensional electron gas 
Eli and recently quantum dots made from isotopically 
purified silicon have shown very long coherence times 
im uni- All this future developments will require more 
and more complex architectures for which the technology 
must be perfectly controlled in terms of reproducibility, 
variability and know-how. In addition, the quality of the 
underlying material must also be addressed as it can be 
an important source of decoherence for quantum devices 
HJ. In that quest, superconducting boron doped sili¬ 
con discovered in 2006 m, can address all of these chal¬ 
lenges at once as silicon technology is by far the most 
mature one for nanoelectronics and for which extremely 
high quality material can be used [ 13 ] ■ Moreover, the 
relatively high normal state resistivity of boron doped 
silicon allows using this material for Microwave Kinetic 
Inductance Detection [HI- 

As a first step towards superconducting silicon based 
quantum engineering, we have fabricated a SQUID (Su¬ 
perconducting Quantum Interference Device) from a sin¬ 
gle layer of superconducting silicon. Our results show a 
flux modulation at low temperature and low field that 
demonstrates the macroscopic nature of the quantum 


wave function of the superconductivity in silicon [TKW5] . 
Our results confirm that boron doped silicon can be ex¬ 
ploited for applied superconducting devices and quantum 
technology. 

The SQUIDs are composed of two weak links (Dayem 
bridges) acting as two Josephson junctions [HI (20] and 
are processed from a single layer of superconducting sil¬ 
icon. The superconducting film has been obtained by 
heavily doping the silicon with boron atoms using the Gas 
Immersion Laser Doping technique [21] . This technique 
consists of shining laser pulses (here an XeCl 308 nm laser 
with a 25 ns pulse duration) at the surface of a silicon 
wafer on top of which molecules of BCI 3 have been pre¬ 
viously chemisorbed. During a laser pulse, the silicon 
melts over a thickness and a time that depend on the 
energy density of the laser. In this melting time, the 
boron atoms diffuse very rapidly into the melted phase 
while the chloride atoms are expelled. At the end of the 
laser pulse, the boron atoms are incorporated in substitu¬ 
tional sites as the Si:B layer is epitaxially grown over sili¬ 
con. Since the number of atoms that can be chemisorbed 
at the surface of silicon is self-limiting, the number of 
boron atoms introduced at each laser pulse is very repro¬ 
ducible and equals 1.2 x 10 14 cm -2 . In order to increase 
the amount of dopants, the overall procedure (gas im¬ 
mersion and laser shot) is repeated (typically 10 to few 
hundred times) to reach active dopant concentrations up 
to 11% at with no boron aggregates [2TI - i?!T| . Such a high 
level of concentration, larger than the solubility limit of 
boron in silicon (~ l%oi), can be achieved thanks to 
the rapid liquid/solid phase transition that quenches the 
boron atoms into the crystalline phase. In order to con¬ 
trol the melting depth, the melting duration is monitored 
through the time resolved silicon reflectivity using a low 
power laser at 675 nm during each laser pulse. This pro¬ 
cedure produces a very sharp doped-undoped interface 
of only a few nm [24] . In our set-up, the surface area 
covered by the laser beam is 2 mm x 2 mm. The study 
of the electronic transport properties down to very low 
temperature has revealed superconducting properties for 
a boron dose larger than 4 x 10 15 cm“ 2 with a maximum 
of the critical temperature below which Si:B thin films 
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show a zero resistance state of 0.7 K for a boron concen¬ 
tration of 10 %at and a thickness of 200 nm [25H27] . 

The superconducting silicon layer we have used for the 
present study was grown with 200 laser pulses and a melt¬ 
ing time duration monitored to 47 ns. With these con¬ 
ditions, we obtained a 80 nm thick layer of boron doped 
silicon Si:B with a boron dose of 2.4 x 10 16 cm -2 , which 
corresponds to a concentration of 5% at. 

The fabrication of the SQUIDs starts with electron 
beam lithography to define the shapes (loops, weak links 
and contact electrodes) of the devices. Then, an alu¬ 
minum hard mask (20 nm) is deposited by lift-off. More 
than 100 nm of the silicon not protected by the aluminum 
layer is etched away by reactive ion etching using fluo¬ 
ride gas and a highly anisotropic recipe. The aluminum 
mask is entirely removed afterwards by wet etching before 
the deposition of Ti/Au contact pads defined by optical 
lithography. Care has been taken to deoxydize the doped 
silicon before the deposition of the contact pads. A SEM 
image of a typical realization is shown in Fig. |T] (left in¬ 
sert). The SQUIDs all have the same loop area and the 
dimensions of the weak links range from 80 nm to 200 nm 
in width and from 100 nm to 500 nm in length. In this 
letter we present the results obtained in a SQUID with 
weak links of 100 nm by 100 nm, but all the devices we 
have measured show magnetic flux modulation. We have 
also checked that the flux modulation in a SQUID loop 
can be reproduced using different silicon batches. 
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FIG. 1. Resistive transition of a superconducting silicon 
SQUID. The SEM image shows the actual SQUID fabricated 
from a bare superconducting Si:B doped layer. It is composed 
of two weak links (Dayem bridge) with dimensions of 100 nm 
x 100 nm for the actual device reported here (see left image). 
Insert : Resistance of the device in the high temperature range 
showing a clear metallic behavior. 

FigjT] shows the resistance of a SQUID as a function 
of temperature. The transition to the superconducting 
zero resistance state occurs at T c = 260 mK as expected 
for such a Si:B layer. The right insert of FigjT] shows the 
resistance of the same device at high temperature up to 
250 A'. The behavior is clearly metallic with a residual 
resistive ratio Rzook/Rwk — 1-4. 

At low temperature, a SQUID critical current I c of 



FIG. 2. Temperature dependence of the measured critical cur¬ 
rent of the silicon SQUID (black square points). The green 
circles correspond to the maximum critical current obtained 
by measuring the magnetic flux modulation and the red tri¬ 
angles result from solving the Ginzburg-Landau equations. 
Left insert : color plot of the phase evolution of the super¬ 
conducting wave function along the device for a phase coher¬ 
ence length of 40 nm and an applied phase fixed to ±3-7t/4 
at both ends of the device. Right insert : IV characteristics 
at 100 mA' showing a strong hysteresis when ramping the ap¬ 
plied DC current. 7 C and I r show the critical current and the 
re-trapping current respectively. 


several /iA (corresponding to a current density of ss 
10 4 A/cm 2 ) is observed in the current-voltage TV char¬ 
acteristics. Below 220 mK the IV response is strongly 
hysteretic (see right insert of Fig|2]) with a very small 
re-trapping current I r . For current larger than the crit¬ 
ical current, the resistance, defined as the slope of the 
TV curve, is 360 Q, which is exactly that of the device 
above the critical temperature. This means that the en¬ 
tire structure turns normal in the dissipative state. This 
behavior is due to a heating effect that propagates from 
the weak links to the rest of the device, as is usually 
the case for such devices. The temperature dependence 
of the critical current measured with no applied field is 
plotted in Fig(2j It shows a rather smooth variation with 
a critical current vanishing at the critical temperature. 

For superconducting Josephson weak links, the am¬ 
plitude of the critical current strongly depends on the 
length of the bridge compared to the superconducting 
Ginzburg-Landau (GL) coherence length £. In the dirty 
limit £ = y/foT, where £o is the BCS (Bardeen-Cooper- 
Schrieffer) coherence length and l the elastic mean free 
path mm- For superconducting silicon at our doping 
level l ss 2 —3 nm and ~ 1000 nm, which gives the zero 
temperature phase coherence length £(0) « 40 — 50 nm 
[ 25 ] . In the actual geometry, the ratio between the length 
of the weak link and the GL coherence length is L/£(0) ~ 
2.5. The critical current of Josephson bridge junctions 
can be obtained by solving the Ginzburg-Landau equa¬ 
tions. Following the code developed by Hasselbach et 
al. [20| . we have numerically solved these equations in 
two dimensions taking into account the exact geometry 
of our device. This calculation is done by assuming a 
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macroscopic superconducting wave function that can be 
written as ip{r) = 00 e ll ^ T 4 where both / and (f> de¬ 
pend on position but ipoc is a constant at a given tem¬ 
perature. The values of /(r) are fixed to / = 1 at the 
left and right sides of the devices (position R and L in 
the SEM image) which means that the amplitude of the 
order parameter is constant and equals to its equilibrium 
value in the electrodes. The simulation is done with a 
certain phase cpu and (f>L imposed at the left and right 
electrodes (usually (j)R = — <^l). We then obtain the am¬ 
plitude and phase everywhere along the device. For the 
simulation, we have assumed a perfectly symmetrical case 
where the critical current in the upper and lower branches 
are equal. The left insert of Fig. [2] shows a color plot of 
the local superconducting phase for a total phase dif¬ 
ference A tf> = 4 >l — 4>R — 37t/2. The simulation shows 
that the main phase drop occurs at the junction bridges. 
Doing the same simulation for various phase differences, 
gives the current-phase relationship of the device and in 
turn, the magnetic field dependence of the critical current 
of the SQUID. The temperature dependence is obtained 
by taking into account the temperature dependence of 
the GL phase coherence length that modifies the ratio 

L/^(T) with £(T) = , The results for the crit- 

yl —T/T c 

ical current at three distinct temperatures is plotted in 
Fig(2] They must be compared, not to the critical cur¬ 
rent obtained at no applied magnetic field, but to the 
maximum of the critical current flux modulation at dif¬ 
ferent temperatures. Indeed, due to remanent and earth 
field contributions, zero applied field does not necessar¬ 
ily mean zero magnetic field. The agreement with ex¬ 
perimental data is very good (see Fig(2j). The general 
temperature dependence obtained here shows a behav¬ 
ior very close to what has been obtained previously for 
similar values of L/£(0) PM5U] , 

One can also estimate the product RnR,wI where Rn 
is the normal state resistance of one bridge junction and 
Ic,wi its critical current. The square resistance of the 
doped silicon film at low temperature can be extracted 
from the normal state resistance of the entire structure 
taking into account the geometry of the SQUID. Since 
the bridge is a square, the normal state resistance of 
one junction is therefore Rn = Rn — 14 D, in agree¬ 
ment with the known resistivity of p ~ 100 /ifl .cm at 
such a doping level and thickness. Considering that the 
SQUID is symmetric, we find RnR,wI — 35^(U with 
Ic,wi = Rl 2 = 2.5 fiA. This value can be compared 
to the superconducting gap A/e ~ 40 [iV that we es¬ 
timate from the critical temperature T c = 0.26 K using 
the BCS relationship A = 1.76fcsT c ■ We then obtain 
RnR,wI — A/e in agreement with existing calculations 

mm 

In Figj3j we show the dependence of the critical current 
with the magnetic field applied perpendicular to the de¬ 
vice, for various temperatures. Each of these curves has 
been obtained by ramping the magnetic field from — 2G 
to +2 G in steps of 0.1 G. The results show a very regular 
oscillation of the critical current as a function of the mag- 
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FIG. 3. Magnetic field dependence of the critical current of 
the DC SQUID at various temperatures from 100 77177 (top) 
to 220777 .K (bottom). The measured modulation (scattered 
points) is compared to the flux dependence obtained from 
the GL equations in the actual geometry for three different 
values of the phase coherence lengths 40 nm, 60 nm and 80 nm 
to account for the temperature dependence (T ~ 100 mK — 
top, 150m .K — middle and 190 mK — bottom) (solid lines). 


netic field with a period of 1 G. This value corresponds 
to a one flux quantum $o = h/ 2e = 2 x 10 -15 Wb in a 
surface area of 20 /xm 2 , which is the surface of the SQUID 
loop 4 /ljto x 5 fj,m. The amplitude is roughly 10% of the 
maximum critical current. This flux (or field) modula¬ 
tion decreases when increasing the temperature and van¬ 
ishes to zero when approaching the critical temperature. 
The solid lines in Fig(3]are the results of the critical cur¬ 
rent modulation obtained by numerically solving the GL 
equations. This has been done for three different values 
of the GL coherence lengths : 40 nm, 60 nm and 80 nm 
corresponding to 100 m/\, 150 m/\ and 190 mK approx- 
imatively. The position of the curve follows rather well 
the temperature dependence (as already seen in Fig{2]) 
and the amplitude of the modulation is well reproduced 
too. 

In conclusion, we have measured the magnetic flux de¬ 
pendence of a DC SQUID device made from a super¬ 
conducting boron doped silicon thin layer. The SQUID 
geometry consists of a 20 fim 2 loop interrupted by two 
nano-bridges forming two Josephson junctions. We ob¬ 
served clear oscillations of the critical current as a func¬ 
tion of the applied magnetic flux. The amplitude of the 
magnetic flux modulation and the temperature depen¬ 
dence of the critical current are well reproduced by nu¬ 
merically solving the Ginzburg-Landau equations. These 
findings are reproducible for samples with various nano¬ 
bridge dimensions and fabricated from different Si:B lay¬ 
ers. These results demonstrate that superconductivity in 
silicon can be described by a macroscopic wave function 
and open routes towards silicon based superconducting 
quantum engineering. The future will be to implement 
more practical devices nans such as superconducting 
resonators and JoFET’s (Josephson Field Effect Tran¬ 
sistors) in which a Josephson current flowing between 
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superconducting drain and source contacts through an 
undoped channel will be modulated by an electrostatic 
gate. Using ultimate technology to reduce the size of the 
channel to only few nanometers, one can expect to ob¬ 
tained tunnel coupling between the drain and the source. 
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